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ABSTRACT 

This paper discusses the issue of using long slender booms as 
pendulous nutation damping device on spinning spacecraft. 
Motivation for this work comes from experience with the 
Galileo spacecraft, whose magnetometer boom also serves as 
passive nutation damper for the spacecraft. Performance 
analysis of a spacecraft system equipped with such a device 
indicates that the nutation time constant of such systems are 
relatively insensitive to changes in the damping constant of 
the device. However, the size and arrangement of such a 
damper raises important questions concerning spacecraft 
stability in general. 


INTRODUCTION 

Most spin stabilized spacecraft are equipped with 
passive nutation damping devices that limit spacecraft 
nutation through on-board energy dissipation. The design of 
these devices is based on well established stability criteria 
for spinning bodies. 1 — ^ When disturbed slightly from its 
position of stable spin, a spacecraft with internal energy 
dissipation will recover faster than one without energy 
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dissipation. This has led to the design of several passive 
devices that are triggered into dissipating energy on board 
of a spacecraft anytime that the spacecraft attitude motion 
is disturbed. Such devices have included simple mass— spring— 
dashpot systems, damped compound pendulum, viscous fluid in 
ring-shaped tubes, etc... The Galileo spacecraft features a 
passive nutation damper that differs markedly from any that 
has been flown to date. As shown in Fig.l, this dual-spin 
spacecraft consists essentially of a rotor, carrying a high 
gain antenna and three long booms, and a stator section that 
houses a probe and carries the scan platform containing most 
of the imaging instruments. The rotor is connected to the 
stator through a spin bearing assembly that allows one degree 
of freedom of relative motion that is controlled by the 
"clock" control loop. The ratio of the rotor spin inertia to 
the vehicle transverse inertia is greater than one (1.4), so 
that the spacecraft spin axis is also its axis of maximum 
inertia. The spacecraft's longest boom - its 
science/magnetometer boom - is also utilized as a passive 
nutation damper by connecting it to the rotor bus through a 
one degree of freedom hinge and a torsional spring and damper 
system, as shown schematically in Fig. 2. Although this 
damper is of the pendulous type, its length (8.6m) far 
exceeds that of any such damper used in past missions. 
Furthermore, the mass center of the boom is outboard of its 
pivot point. This, again, is unusual for pendulum dampers. 

In the remaining part of this paper, the effectiveness 
of this design is explored and compared with that of 
traditional dampers. This is done essentially by examining 
the shape of the damper time constant versus damping constant 
curve as well as by studying the overall attitude dynamics of 
the spacecraft in the presence of such a large boom. 
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EQUA TIONS OF ATTITUDE MOTION 


To derive the dynamical equations of attitude motion of 
the spacecraft, the model shown in Fig. 3 is used. The dual- 
spin nature of the spacecraft is ignored, and the system is 
assumed to be made up of a main rigid body A and a boom B. A* 
and B are mass centers of A and B respectively, and S* is 
the mass center of the combined system. The following 
simplifications are also made: 

* S* lies on the spin axis Z of the spacecraft and 
remains fixed in body A at all times; 

A , B and S all lie on a plane containing the Z 
axis; 

* when P = 0 , the central principal axes of the 
system for S are parallel to a^,a 2 , and 33 ; 
bi,b 2 , b 3 are parallel to central principal axes of 
B for B*, and, for body B, = i 3 , while I 2 = 0, 
where I indicates moment of inertia. 

The equations of attitude motion of this system, as derived 
using AUTOLEV 6 ' 7 , are: 

§(Li}Sinf} + zo ) + (LqcosP + yo)\m B + (y][ + zfyniA + (/f + + 

{[(LffSinP + z Q )L B sinP + {LgcosP + yo)L B cosp)m B + I D )ii 4 - 
{[(«! + u 4 )L B sinP + zomi]«i + 

(L B u 2 sinP - L b u 3 cosP + z 0 U 2 -you 3 )u 2 + 0<i + u 4 )L B sinPu 4 )(L B cosp + y 0 )m B + 

{-[(?<i + u 4 )L b cosP + }’o/q]»i + \L B (u 2 sinp - u 3 cosp) + zqii 2 - yo/oj'O - 
(ui + u 4 )L b cosPu 4 } ( L B smP + zo)m B - [[(zaU 2 - yAii 3 )u 3 - yAufizA + 

[0/U<2 - yAU 3 )u 2 + ZAH\}yA}tnA + (u 2 cosp + u 3 sinp)(u 3 cosP - u 2 sinp)I B - 

( l 2 u 2 + 7 23 ,< 3)«3 + Un u 2 + ^ ll i) u 2 = 0 
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\(L B sinP + zo) m B + (l 2 + l D sin 2 p) + m A z^\ii 2 - 
[{Lficosft + yo)(L B sinp + z B )m B + l B sinfkosf$-Izs + fn A y A z A ] j '< 3 + 
{[(/<i + u^Lncosfi + >’oMi]« 2 + [o<i + ii4)L B sinP + zqHi]/< 3 + 

( u 2 cosP + ujsinP )L B it4 ) (Lbs inf} + zo)m B + ([(uj + 114)112 — ii\U2]sinP + 
[(ill + 1(4)112 - U\U2 ]cosP)I D sinP + ( u 2 cosP + u 2 sinP)(ui + U4)I°sinp- 
(f 23 l< 2 + / 3 /<3)»l + 0/1MJM2 + ZAll\U3)m A Z A + /f«if< 3 = 0 


[(LrcosP + yo)(L/)SinP + zo)mo + I B sinPcosP -1 22 + m A y A z A ] u 2 - 
[( LqcosP + jo) nig + 1 B cos 2 P +/ 3 + m A y\ ];) 3 + 

{[(«i + U 4 )LbcosP + >w<i]»2 + [(«i + i( 4 )Li)sinP + z 0 /<i ]/< 3 - 
L B cospu 2 U4 + L B sinPu2tt4)(L B cosP +y 0 )m B + {[(//j + 114)113 -uiii 3 ]cospsinp + 
+ ii4)u 2 — U\ii 2 ]cos 2 P)I B + (u 2 cosP + ii2Sinp)(ii\ + 114)!^ cosP — 

(/23U2 + / 2 " 2)"1 + (yAllll(2 + ZAU\ll3)m/iy A + rfinu 2 = 0 


{[(L B sinP + z 0 )L D sinp + ( L B cosP + y 0 )L B cosp]m B + I D )ii x + (/ D + m B L 2 D )ii4 - 
f[((«i + U4)L B sinp + z 0 hi)ui + (L B sinpu 2 - L B cospU2 + m2 - yoit 3 )u 2 + 

(ill + ii4)L B sinPit4\L B cosP - [((//] + U4)L B cosP + yQUiju\ - 
(~L B ii2CosP + L B u 2 sinP + zqii 2 - y 0113)113 + 

(u\ + ii 4 )L B cosPii 4 \L B sinP}m B + (u 2 cosp + U2sinp)(ii2CosP- u 2 sinp)I B + 

kp + 01/4 = 0 


where uj(l - 1,2,3) are the components of the angular 
velocity of A along a^, u^= P, m represents mass, I 
represents moment of inertia, k is spring stiffness, c is 
damping constant, and the dimensions y 0 , z 0 , y A , z A , L B are as 
shown in Fig. 3. 



NUT AT TON DAM PER TIME C ONSTANT 


A known equilibrium point of the system corresponds to 
the pure spin condition. That is the solution 

Ul = u 2 = 0 , U 3 = Q (const.), U 4 = 0 , and (i = 0 . 

This solution does satisfy Eqs . (1-4) provided that 

1 2 3 = m/i z /\Ta + /riflzo (Lb + To) ( 5 ) 

a condition that is indeed satisfied by the inertia related 
simplifying assumptions given earlier. When the full 
nonlinear dynamical equations given as Eqs (1—4) are 
linearized about the pure spin solution, the result is a set 
of first order differential equations that has the form 

Bx T = Ax T ( 6 ) 


where 


x = [u\ ui 1/4 p] C7) 

and A and B are 4 by 4 matrices with the following elements: 

A 1 1 - /l 13 = An ~ An — A14 - A3 1 = A4 1 = A 41 = A 44 = 0 (8) 

An - + To ) 2 + ma(T 2 - z a) + I D - i} + ^3]^ 

A 14 = ~[l D + + jo) }& 0 ^) 

A 2 1 = -Ilf - 7 3 m A z A + m B z 0 Cl 1) 

A 32 = ~[l D + ui[)Lb(Lb + To)]^ (12) 

A 3 3 = — (T, A 43 = 1 ( 13 ) 
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(14) 


A 34= -j m B L B {L B + >’o)^ 2 + Jfc] 

B] \ = [l\ + 1 D + m A (yl + i\) + m B [(L B + y 0 ) 2 + 2 fl) 

B 12 = #14 = #21 = B 2 3 = B 2 4 = Z?32 = Z ? 34 = fl 4l = Z ? 42 = Z? 43 = 0 
B\3=[l B + m B L B {L B + jo)] 
b 22 = Vi + rn A z\ + 

/? 3 i = [f B + m B L B (L B + jo)] 

#33 = f B + m B L\ 


( 15 ) 

( 16 ) 

( 17 ) 

( 18 ) 

( 19 ) 

( 20 ) 


#44=1 (21) 

The eigenvalues of the matrices A and B are found to have 
negative real parts for inertia property values corresponding 
to all mission phases of the spacecraft. Hence, the pure 
spin solution is a stable solution. The nutation angle time 

constant is the negative reciprocal of the eigenvalue 
corresponding to uj_ or U 2 . 



Fig. 4 shows Galileo's nutation angle time constant 
plotted against the damping constant of the passive nutation 
damping device on board. The case shown corresponds to a 
damper spring stiffness of 335 N. m/rad, and spacecraft 
inertia property values near the beginning of the mission. 

Two important facts emerge from this plot. First, there 
is an optimum damping constant corresponding to a given 
spring stiffness. The most remarkable thing about the curve 
shown is the fact that it is so flat; especially near the 
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minimum time constant value. This means that there is a wide 
range of values of a for which the time constant changes very 

little. This result is in great contrast with what obtains 
for traditional passive dampers, where such plots are not 
flat at all, and "tuning" of the damper is almost always a 
necessity if one desires reasonably small time constants. 

This relative insensitivity of the time constant to O is 

particularly appropriate for interplanetary missions. This 
is because the viscosity of damper fluids is generally very 
sensitive to temperature, and, therefore, the damping 
constant can be expected to vary widely during a long 
interplanetary flight that takes a spacecraft through varied 
environments. It is thus advantageous to have a damper, 
whose performance will not be degraded by the inevitable 
fluctuations in damper fluid viscosity. 


As exemplified by the design and performance analysis of the 
Galileo passive damper system, the use of long booms as 
nutation damper for spin stabilized spacecraft introduces a 
new and important advantage over traditional damping devices. 
It renders the system nutation angle time constant 
practically insensitive to the device damping constant, 
thereby drastically reducing the need for "tuning" of such 
dampers , The mam disadvantage of such a large device is 
that it becomes an important factor in spacecraft stability. 
Furthermore, because of the small relative damper 
displacements that are to be expected from this design, 
factors such as stiction become important in the evaluation 
of the damper's performance, and may impose thresholds on the 
amount of nutation that can be damped out. 
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